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We calculate general moments of the lepton energy spectrum in inclusive semileptonic B — > 
X c tv decay. Moments which allow the determination of m\ s and Ai with theoretical uncertainties 
Aml s ~ 0.04GeV and AAi ~ 0.05GeV 2 are presented. The short distance 15 mass is used to 
extract a mass parameter free of renormalon ambiguities. Moments which are insensitive to mb and 
Ai and therefore test the size of the 1/mf matrix elements and the validity of the OPE are also 
presented. Finally, we give an expression for the total branching ratio with a lower cut on the lepton 
energy, which allows one to eliminate a source of model dependence in current determinations of 
\V c b\ from B — > X c lP decay. 



Inclusive semileptonic B decays have been examined 
in recent years using an operator product expansion 
(OPE) || . This OPE coincides with the parton model 
in the limit m& — > oo with nonperturbative corrections 
suppressed by powers of Aqcd /mb- There are no nonper- 
turbative corrections at order Aqcd/w&, and the order 
(AQCD/wb) 2 corrections are parameterized by two ma- 
trix elements 10. 
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Here h v is the heavy quark field in the effective theory. 
Physically these parameters correspond to matrix ele- 
ments of the kinetic and chromomagnetic operators in 
the heavy quark effective theory [Q. To second order in 
Aqcd/to6 the semileptonic decay rate for B — > X c £p in 
the OPE is 
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g{p) = -9 + 24p-72p 2 + 72p 3 -15p 4 -36 y o 2 log /9 . 

This inclusive decay can be used to measure the magni- 
tude of the CKM matrix element \V c b\. The results of 
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this analysis yields a value in good agreement with de- 
terminations from exclusive decays S 

\Vcb\ = (40.4 ± 0.9 exp ± 0.5 AljA ± 0.8th) x lO" 3 . (4) 

The first uncertainty is experimental, the second is due to 
uncertainties in the values of mb and Ai, while the third 
is from unknown higher order terms in the OPE and 
unknown perturbative corrections. One of the system- 
atic uncertainties is due to modeling the spectrum below 
the required cut on the charged lepton energy spectrum 
and therefore is also theoretical in nature. For machines 
running on the T(AS) resonance, this experimentally re- 
quired cut is typically at 0.6 GeV. 

The measurement of \V c b\ from inclusive decays re- 
quires that these decays be adequately described by the 
OPE formalism. In particular, the possibility of devia- 
tions from the OPE predictions due to quark-hadron du- 
ality have been raised ||. However, in order to compare 
the OPE predictions with data, one also has to define 
how uncertainties from 1/m^ corrections are estimated. 
These uncertainties are hard to quantify reliably, as the 
only information on the size of matrix elements of the di- 
mension six operators comes from dimensional analysis. 
The important question to answer for a reliable deter- 
mination of \Vcb\ from inclusive decays is therefore how 
well the OPE fits the data. Thus, in this paper we define 
duality violation to be the difference between the OPE 
predictions using dimensional scaling of the 1 jm\ matrix 
elements and the experimental data. 

It is the purpose of this paper to demonstrate how to 
reduce the theoretical error on the extraction of |T4f,|, m\, 
and Ai from semileptonic B — ► X c tv decay. We elimi- 
nate the uncertainty due to modeling the spectrum below 
0.6 GeV by calculating the dependence on this cut. We 
demonstrate how to improve the precision of the determi- 
nation of the parameters mb and Ai using measurements 
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of generalized moments of the lepton energy spectrum. 
We also calculate several moments which are insensitive 
on the actual values of m& and Ai and therefore allow 
one to test the validity of the OPE and therefore our 
definition of duality violation. 

Past attempts to measure mt, and Ai from the lepton 
energy spectrum of B — > X r tv used the moments pre- 
sented in 
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where £ max = (m B — m? D ) /2m B . These moments were 
calculated to third order in the Aqcd /mb expansion 
H and to order a 2 [3 || in perturbation theory. The 
pole masses were expressed in terms of the heavy meson 
masses and the parameter A via the relation 
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It is well known that the pole mass and the parameter 
A have a renormalon ambiguity p0| . In every physical 
observable, such as these moments, this ambiguity is can- 
celed by corresponding terms in the perturbative series 
p[ . This means, however, that a value for the pole mass 
(or A) only has a meaning if specified at a certain order 
in perturbation theory. While this is not a problem in 
principle, it is much more convenient to work with a short 
distance mass which does not exhibit such a renormalon 
ambiguity. Using a short distance mass instead of the 
pole mass also tends to give a perturbative expansion 
which converges more rapidly. 

In this paper we calculate the general moments of the 
lepton energy spectrum 1 
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to order (Aqcd/to6) 3 and a 2 (3$ using the short distance 
IS mass and systematically examine what information 
can be obtained from these moments. In particular, we 
identify sets of moments which allow mb and Ai to be 
measured with uncertainties considerably smaller than 
Ri and i?2- We also find moments which are insensitive 
to the values of toe, and Ai and therefore directly test 
the reliability of the theoretical error assigned to the ex- 
traction of TOt, and Ai and the theoretical treatment of 
inclusive B decay using an OPE. 

The moments are calculated from the lepton energy 
spectrum, which has been calculated to order (A/mt,) 3 



1 In general, non-integer moments can lead to new cuts in the com- 
plex plane which have to be dealt with carefully. This, however, 
is not a problem for this paper. 



[U, and to order a 2 (3 perturbatively || The cal- 

culated moments depend on the matrix elements of the 
higher dimensional operators and the pole masses of the 
b and c quarks. To avoid the renormalon ambiguity, we 
rewrite the expressions for the moments in terms of the IS 
mass, which is related to the b quark pole mass through 
the relation Fl, p| 
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where i — log[///(m&a s CF)] and Cf = 4/3. The depen- 
dence on the pole mass of the charm quark is eliminated 
through the relation 
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where p\ and are matrix elements of local dimension 
six operators ||. In this relation we use the fact that 
ifiB — m\ s ~ Aqcd by expanding in the parameter A\s 



Ais = m B - m\ s . 



(10) 



It is important to note that in contrast to A, which re- 
lates the pole mass to the meson masses, the parameter 
Ais is not an infrared sensitive quantity. Using m\ s ne- 
cessitates introducing a modified perturbative expansion 
in order to ensure the cancellation of renormalon ambi- 
guities When calculating in the IS mass scheme 
the order a s n perturbative corrections coming from the 
mass transformation Eq. (|^) are counted using the pa- 
rameter e n_1 while a™ corrections in the decay rate are 
counted as e n . The parameter e determines the order in 
the modified perturbative expansion. 

The lepton energy spectrum has been measured by 
the CLEO collaboration using single and double tagged 
charged lepton data samples |15|, [16| . The double tagged 
data uses the charge correlation between primary and 
secondary leptons to eliminate the background from sec- 
ondary leptons. While this type of analysis gives the 
smallest experimental uncertainties, it does require a cor- 
rect identification of the primary lepton. For this reason 
a cut on the lepton energy at Ee — 1.4 GeV was employed 
by CLEO jl5| . The double tagged data sample does not 
require such a cut, and one only needs a lower cut on 
the lepton energy at Eg = 0.6 GeV to eliminate fake 
pion signals. However, this increases the uncertainties in 
the moments by about a factor of three Jll| . With the 
large data samples available to the B factories it should 
be possible to use single tagged data samples to measure 
the moments with the required precision. For each search 
of optimal moments we therefore present results relevant 
for the two analysis techniques by restricting the lepton 
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energy to lie in the conservative regions above 800 MeV 
and above 1.5 GeV, respectively. 

To compare with experimental data at the required 
precision, we also calculate the corrections due to the 
electroweak electron radiative corrections |l7|, the lead- 
ing order corrections proportional to IKib/VchI 2 , and the 
corrections due to a boost from the rest frame of the B 
meson to the lab frame. For the boost corrections we as- 
sume a mono-energetic B meson of energy E — My /2 
appropriate for the CLEO analysis. Boost corrections for 
the asymmetric B factories are not considered here and 
should be taken into account in the experimental anal- 



ysis. In this case, the boost corrections calculated here 
should not be included. 

Before we present the results of the general moments 
search, we can use our result to obtain the total branch- 
ing ratio with a cut on the lepton energy. As explained 
earlier, such a cut is required experimentally to avoid 
fake lepton backgrounds. Here we present the branch- 
ing ratio as a function of the lower cut on Eg. Since the 
dependence on the cut is quite complicated, we will not 
present an analytic equation, but rather give an approx- 
imate formula, which deviates from the exact result by 
less than 1% for < E cut < 0.8 GeV 
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Here E cut = E cut /m B and we have used a s = 0.22, 
A, = 25/3 and A 2 = 0.12 GeV 2 . The dependence of 
the uncertainties in this expression from perturbative and 
1 jm\ corrections is negligible over the range of E cut con- 
sidered, and can thus be taken as constant. 



For the measurement of mj and Ai, we compare the 
resulting uncertainties on the extracted parameters with 
those obtained using the moments Ri and R2 defined 
in (||). To facilitate this comparison, we present these 
moments in terms of the IS mass 
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R 2 = #[0,1.7,0,1.5] = 0.6578 
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The last two brackets in each moment are the electroweak CLEO collaboration 
and boost correction, respectively. Using the central val- 
ues and statistical correlation matrix V measured by the #1 = 1 7831 #2 = 6159 
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together with \V ub /V cb \ = 0.1 ± 0.03 and a s = 0.22, we 
obtain 

Ais = [0.47 ± 0.10 exp ± 0.02 Vui ± 0.02 e ± 0.07 m3 ] GeV 
Ax = [-0.16 ± 0.1l exp ± 0.02 Kt ± 0.02 e ± 0.07 m3 ] GeV 2 , 

where the perturbative uncertainties are estimated by 
varying m b /2 < /i < 2m b and using the O(e 2 f3o) term as 
a crosscheck. The uncertainty due to the unknown size of 
matrix elements of 1 jm\ operators is estimated by using 
the method of |8[ [ll| to ensure a valid comparison of our 
suggested observables with Ri and R^'- We impose the 
relation due to the vector, pseuodoscalar mass splitting 
at third order § 

92 - T2 - r 4 = 

Km 2 B Am B (mp + Ais) - m 2 D Amp(m B + Ais) 
m,B + Ais - k {mo + Ais) 

and use the result from the vacuum saturation approx- 
imation which predicts a positive value of p\ fH| . The 
unknown matrix elements are then randomly varied be- 
tween (±500 MeV) 3 . The random values are drawn from 
a flat distribution, since there is no known preferred value 
of the individual matrix elements. 

We would like to make a few comments about the ex- 
pressions for R\ (|l^ ) and i?2 (|l3|). First, the expressions 
differ from the ones presented in [Q in their perturbative 
terms and the terms containing Ai. This is due to the 
use of a short distance mass in the present paper, as sup- 
posed to the pole mas used in Q. Furthermore, in both 
moments the O(e 2 (3o) term in the perturbative series is 



larger than the order e term. This, however, is due to a 
cancellation of the O(e) terms in the ratios considered, 
and does not indicate a poorly behaved perturbation se- 
ries. In both the numerator and the denominator of these 
expressions the perturbation series is well behaved. 

Translating the obtained value of Ais into a value of 
the b quark mass we find m\° = 4.84±0.13, adding the er- 
rors in quadrature. This result is in agreement with both 
sum rule and other moment analysis extractions |^0j , with 
comparable uncertainties. Using the definition of the 
general moments given in (Q), one can search for com- 
binations of moments which allow for a determination 
of m\ s and Ai with smaller uncertainties. For certain 
values of n, m, Eg 1 and E( 2 , the OPE results can not 
be trusted any more since the convergence of the OPE 
breaks down. Using dimensional analysis to obtain the 
size of the matrix elements and requiring that the 1 / m 2 
and l/f7i| contributions do not exceed 10% and 3% of 
the leading order result, respectively, we restrict ourself 
to the parameter space 

m< 3 , n < 3 , 0.8 GeV < Eff 1 < 1.8 GeV (14) 

in our search to ensure a well behaved OPE. 

By simultaneously minimizing the size of third order 
corrections and maximizing the linear independence of 
the moments we find sets of observables R a and R b which 
minimize the theoretical uncertainty in extracting Ais 
and Ai. We do not try to minimize the size of the per- 
turbative corrections. Since the full two loop correction 
to the lepton spectrum is not available, it would be mean- 
ingless to minimize the size if the e 2 /3o piece. Consider 
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For this moment, the only available data is CLEO's dou- 
ble tagged lepton spectrum, and the uncertainties are 
too large to be useful for this analysis. Thus, to il- 
lustrate how well the parameters A iS and Ai can be 



extracted from these moments we use the hypotheti- 
cal data R { a ] = 1.0082 and R^ = 0.9266. Using 
again \V ub /V c b\ 2 = 0.1 ± 0.03 and a s — 0.22, this leads 
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FIG. 1: Comparison of the error ellipses due to the 1/mf 
uncertainties for the various sets of moments. Only the central 
value of the black, solid ellipse has meaning, for all the other 
ellipses only their relative size is important. 

to Ais = [0.59 ± 0.02 Vub ± 0.03 e ± 0.01 m3 ] GeV, and 
Ai = [-0.16 ± QM Vub ± 0.02 e ± 0.03 m3 ] GeV 2 . Adding 
all errors in quadrature, this leads to a theoretical error 
on m\ s of ±40 MeV, with the largest error from the per- 



turbative uncertainties. To further reduce the error, a 
full two loop calculation for the lepton energy spectrum 
would be required. The largest uncertainty on Ai is due 
to the error in \V u b\, for which we have assumed a very 
conservative error. Future measurements of |V^b| should 
lower this uncertainty considerably J2l]j . In Fig. |l| we 
compare the resulting 68% confidence level ellipses in the 
A15 — Ai plane with the one obtained from R± and i?2- To 
estimate how the experimental error from these new mo- 
ments compares with that from the previous extraction, 
we scale the measured correlation matrix for R\ and R2 
by the central values to keep the percentage error equal. 
This leads to slightly increased experimental uncertain- 
ties compared with those using R\ and i?2- Thus, to 
minimize the overall uncertainty, the new moments have 
to be measured with smaller uncertainties than R\ and 
i?2 jl(|. Considering that this measurement was based 
on only 2/6 -1 of data, this seems feasible. 

Next we turn our attention to moments that are insen- 
sitive to the values of nib and Ai and are therefore well 
suited to test the underlying assumptions of the OPE 
and thus for duality violations. Consider 
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Using the same input parameters as before, and very con- 
servative errors for the two HQET parameters A15 = 
0.6 ± 0.2, Ai = —0.3 ± 0.3, this moments is predicted to 
be 

L>i = 0.5400 [l ± 0.0010 Al ± 0.0003 Ais ± 0.0006y„ b 
±0.0002 £ ± 0.0009 m3 ] . 

Thus, even with little information on the size of A15 and 
Ai, the numerical value of this moment is predicted to 
better than 1% in the OPE, given that the method to 
estimate the 1/m^ uncertainties is reliable. A measure- 
ment of this moment can therefore directly test these un- 
derlying assumptions, and any deviation from the OPE, 
even if below the % level, should be detectable. In the ap- 
pendix we give several other moments which have similar 
theoretical uncertainties. While each of these moments 
individually tests the accuracy of the OPE predictions, 
it is best to measure several of these moments. This will 
give important information on the validity of using the 
OPE for other inclusive measurements, such as the ex- 
traction of |V c fc|- 



We have presented observables that allow one to im- 
prove on the measurement of the parameters ml s and 
Ai. The theoretical error in the extraction of both ob- 
servables is improved, with the largest remaining uncer- 
tainty being from the perturbative expansion. Thus, to 
further reduce the error, a full two loop calculation of the 
lepton energy spectrum would be required. We have fur- 
ther presented moments which are insensitive to the size 
of these two parameters and therefore allow to directly 
test the underlying assumptions of the OPE. Both mea- 
surements are important to reduce and gain confidence 
in the theoretical uncertainties present in the inclusive 
determination of \V c b\- The method of searching a gen- 
eral set of moments could also be extended to different 
decay distributions, such as the hadronic invariant mass 
distribution. 
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APPENDIX A: ADDITIONAL MOMENTS 

In this appendix we present several additional moments which either minimize the theoretical uncertainties on the 
nonperturbative parameters Ais and Ai or which are insensitive to these parameters. For each of these categories we 
give moments subject to the constraint Eg > 1.5 GeV and moments in which we consider the full range of Eg. 



1. Moments to extract Ais and Ai with no restrictions on Ei 
Moments R[1A, 1.3, 0.8, 0.9] and #[1.6, 1.4, 0.9, 0.8] 
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Using the hypothetical data R y a ' = 0.9096 and R\ ' = 1.5666 we find 



Ais = (0.62 ± 0.01y ub ± 0.04 e ± 0.01 m s) GeV , Ai = (-0.19 ± 0M Vub ± 0.06 e ± 0.03 m3 ) GeV 2 . 



(A3) 



2. Moments to extract Ais and Ai, restricted to E e > 1.5 GeV 
Moments #[0.7,1.7,2,1.5] and #[0.9,1.6,0,1.7] 
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T4 



m 



m 



v, 



cb 



0.22 + 3.7 



Ais 
m B 



0.0076 + 0.014 



B 

Ais 

m B 



e 0.0003- 0.001- 



A 



is 



m B 



0.0016- 0.006 



Ai< 



m B 



#[ 3) = #[0.9,1.6,0,1.7] = 2.2041 



1 + 0.14^ + 0.554¥ + 0.40^- + 2.58^ + 1.7^# + 4.9 AlsAl 



m B 



mi 



rn 



+15.6^-2.74 

m B m B 



2.94_ + l.o4 

m B m B 



0.9- 



r 2 



2.2- 



B 
T3 



2.6 



mi 



mi 



mi 



m 



B 



mi 



e 0.0001 - 0.001- 



A 



is 



m B 



-e| LM 0.0006 



Vub 



v, 



cb 



0.11 + 2.3^ I + (0.0044 + 0.01o4^| + (0.0013-0.0024^ 

m B J V rn B J \ m B 



Using the hypothetical data #i 3) = 0.2955 and #<; 3) = 2.2908 we find 

Ais = [0.64 ± 0M Vub ± 0.01 £ ± 0.02 m3 ] GeV , Ai = [-0.18 ± 0.05y ub ± 0.04 e ± 0.02 m3 ] GeV 2 



(A4) 



(A5) 



(A6) 
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Moments R[0.8, 1.6, 0, 1.7] and #[2.5,1.6,2.9,1.5] 



#W = #[0.8, 1.6, 0,1.7] = 2.0712 



1 + o.l6^1£ + o.574¥ + 0.55^ + 2.79^ + l.l^f + 5.4- A f Al 

™b 

Ais 



nig m?, ' " J "rn 2 R ' ' ' Tod to 3 



(A7) 



to. 



.l4V - 2.7^2- + 1.24V - 0.64V + 2.34V + 2.84V - e f 0.0001 - 0.001 

V m B 



-e| LM 0.0008 



to. 



Kb 



to. 



TO, 



TO, 



TO i 



TO, 



K 



cb 



i?[ 4) = 



f 0.19 + 1.7^ | + ( 0.0046 + 0.01o4^ ) + (0.0008-0.0034^ 
V m B J \ m B J \ m B 

.Ai.<? A? c „ Ai , „_ A 2 „ n A 3 s 0/ ,Ai S Ai 



i?[2.5,1.6,2.9,1.5] = 0.6933 [l- 0.09^ -0.30^ -0.40^ -1.65^§-- 0.8^ -3.6^-^ (A8) 



TO , 



Hi, 



TO, 



TO, 



-8.7^ + .9^V + 1-4^- - 0.8^L + 0-2^ - 1.44- - 1-64- + * ( 0.0002 + 0.0002^ 

TOd TOg TOp TO^; TOp TOp TO^ \ TTl B 



B 



B 



B 



"B 



B 



+e| LM 0.0005 + 



Kb 



0.12 + 1.3^ ) - (0.0031 + 0.006^ ) - (0.0003-0.0034^ 

m B J V m B J \ m B 



Kb 

Using the hypothetical data #i 4) = 2.1558 and #[ 4) = 0.6788 we find 

Ais = (0.64 ± 0.01y u6 ± 0.01 £ ± ±0.02 m3 ) GeV , Ai = (-0.19 ± 0.06y ub ± 0.02 £ ± ±0.02 m3 ) GeV 2 . (A9) 



3. Moments with small dependence on Ais and Ai with no restrictions on Et 
Moment #[0.8,1,0.1,1.3] 

D 2 = #[0.8,1,0.1,1.3]= 1.7587 [l + 0.0l4^ + 0.1l4^ + 0.005^ + 0.97^|- + 0.34¥ ; + l-2^V^- (A10) 



m B 



TO, 



TO, 



TO, 



TO, 



TO 



B 



+ 4. ^_l. iL 



1.74V + 0.34V - l.l4V + 1.04V + 1-04V - e (0.0010 + 0.0024^ 

m B m B TO B m B m B V TO B 



"B 



B 



B 



-e^ LM 0.0004 + 



Kb 



Kb 



0.02 - 1.2^ | + (0.0055 + 0.0084^ ) + (0.0015 + 0.0004^ 
m B J V m B ) \ m B 



This leads to 



D 2 = 1.7784 [1 ± 0.0014 Al ± 0.0016 Als ± 0.0007y u() ± 0.0004 e ± 0.0009 m3 ] 



(All) 



4. Moments with small dependence on Ais and Ai, restricted to Eg, > 1.5 GeV 
Moment #[0.7,1.6,1.5,1.5] 



D 3 = #[0.7,1.6,1.5,1.5]= 0.5256 



1 _ . 04 Al£ - 0.24% + 0.09^ - l.ie^- - 0-7% - 2.0^i (A12) 
mg rot m R m B Tod to„ 



-7.8^1 + 2.6-^ + 2.2^- - 0.2^ + 1.34 - 1.04 - 1.24 + e ( 0.0001 + 0.0001^ 

TOo TOr> TOd TOd TOd TOd TOd V TOfi 



+e| LM 0.00005 



Kb 



Kb 



0.13 - 3.o4^ I - ( 0.0030 + 0.007^ | - ( 0.0018 - 0.0008^ 

m B J \ m B J \ m B 



This leads to 



D 3 = 0.5166 [1 ± 0.0014 Al ± 0.0042 Ais ± 0.0013^ ± 0.00005 £ ± 0.0012 r , 



(A13) 



8 



Moment #[2.3,1.6,2.9,1.5] 



D 4 = R[2.3, 1.6, 2.9, 1.5] = 0.6108 

1.8 



.7.9^ + 2.5^ 

m B m B 



P2 



-BLM 



0.0002. 



v m 



V, 



cb 



0.05 — - 


A 2 

- 0.26^r 


m B 




-0.3^- 

m% 


m B 


3.5^) 


- ^0.0027 







0.02 



Ai 



1.17- 



0.8 



4k_ 2 .6^i (A14) 



1.0- 



T3 



1.2 



T4 



B '"B 

37^ V 

TO B J 



e 0.0001 + 0.0003 



Aig 



Al£ 



0.0013 - 0.003 



This leads to 



L> 4 = 0.6016 [1 ± 0.0034 Al ± 0.0042 Ais ± 0.0017y ub ± 0.0002 e ± 0.0011 m3 ] 

I 



(A15) 
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